Mixing of two-electron spin states in a semiconductor quantum dot 
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We show that the low lying spin states of two electrons in a semiconductor quantum dot can be 
strongly mixed by electron-electron asymmetric exchange. This mixing is generated by the coupling 
of electron spin to its orbital motion and to the relative orbital motion of the two electrons. The 
asymmetric exchange can be as large as 50% of the isotropic exchange, even for cylindrical quantum 
dots. The resulting spin mixing contributes to understanding spin dynamics in quantum dots, 
including light polarization reversal. 



An electron spin in a semiconductor quantum dot (QD) 
is an attractive qubit for quantum computing pj: the 
spin in the ground orbital state can have long coherence 
time 0; a single qubit can be initialized or read opti- 
cally by transient electron-hole pair excitation giving a 
negative trion X~ 0, E|; and the manipulation of the 
spin exchange between neighboring spins can be the ba- 
sis for two-qubit gates . A detailed picture of correla- 
tions between spins in QDs is essential for understanding 
the spin dynamics. The dominant interaction between 
two electrons (e-e) is the Heisenberg-like spin-symmetric 
Js\-S2 (symmetric exchange), which conserves the total 
spin S=S\+S2- Additional spin-asymmetric e-e interac- 
tions (asymmetric exchange) do not conserve S and thus 
decrease the fidelity of gate operations. 

Among a number of recent experiments giving informa- 
tion about spin dynamics are those involving an optical 
polarization reversal 0, 0] . For them, it has been sug- 
gested that this effect results from spin flipping due to 
electron- hole (e-h) exchange in QDs with lateral asym- 
metry. However, these experiments require strong spin 
mixing, inconsistent with e-h exchange alone Q . 

Spin-orbit (s-o) interactions play a key role in under- 
standing mixing of spin states. They arise from effective 
magnetic fields created by the orbital motion of electrons 
p| . Electrons in QD ground states with dominant s com- 
ponents have small orbital angular momentum and thus 
small s-o coupling. A number of experiments of interest 
involve excited electrons in excited states of the QD. Lin- 
ear combinations of nearly degenerate excited states in 
a plane (e.g. p x - and p y -like) can give rise to 2D orbital 
motion with an effective magnetic field perpendicular to 
the plane, and thus to large s-o coupling. This is analo- 
gous to the L-S coupling in atoms [6|. Thus, relatively 
symmetric QDs (e.g. cylindrical) can have significant s-o 
effects, as we show here. 

There are three sources of s-o coupling that lead to the 
mixing of spin states. The largest two contributions arise 
from the kp mixing between the conduction and valence 
bands near the zone center, as described in the effective 
mass approach Q. We derive them by treating the po- 
tentials from the structure and from the e-e Coulomb 
repulsion on the same footing with kp terms, using the 
Kane model ; $j . We have in mind QDs with a strong con- 
finement to a single state £(z) along the growth axis e z , 



and a weaker confinement in the transverse directions, 
which give the electron states (fri(r)=£(z)<pi(p). 

The first contribution to the s-o coupling, h v , arises 
from the structure potential V(r) of the QD. It gives a 
single-electron s-o coupling of the form 0: 
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where p z is not present due to the strong vertical confine- 
ment (for a single state £(z), (£|p z |£)=0). The first term 
in Eq.Q is the usual Rashba coupling j v (e z xp- 1 ), where 
7 V =^Y (£.\d r y\£) , associated to asymmetry in the growth 
direction [lCf. The second term is important for excited 
states whose main components are inversion-asymmetric 
(p-like), where it gives the dominant s-o coupling, inde- 
pendent of structure or bulk inversion asymmetries. This 
term vanishes in the QD ground state, whose main com- 
ponent is inversion-symmetric (s-like). 

The second contribution, h c , arises from the interac- 
tion of each spin with the orbital motion of the other. 
We have obtained it within a two-particle k-p approach 
for electrons interacting through the Coulomb potential 
Uc( r r) [TTlfT^. For k =1,2 and r r =r\— r%, we have: 



s k = (-l) k "/ s CV rr Uc x p k ) • Sfe 



(2) 



The coupling h v from Eq. (JIJ is analogous to the Pauli 
s-o interaction, while h c from Eq.J2J is analogous to the 
Breit-Pauli spin-relative orbit coupling The Pauli 
and Breit-Pauli couplings in vacuum or in atoms are rel- 
ativistically small, due to the large energy gap 2moC 2 be- 
tween electron and positron bands, whereas the present 
gap Eg is smaller, giving larger s-o couplings. 

There is also a smaller contribution, h B , from the Dres- 
selhaus coupling due to the lack of bulk inversion symme- 
try |13| . It arises from the mixing of the conduction band 
with the remote upper bands and it gives a single-particle 
s-o coupling in the form h B -s—j^e al35 p a (p'^—pg)s a , 
with indices denoting crystal symmetry axis. In the 
QDs with strong vertical confinement consider here 
(p z2 )^>(p 1 ~' 2 ) , thus the effective Dresselhaus coupling con- 
tains only the transverse components h B ^^^ B (p x , —p y ), 

with 7 S =7fM0. 

We use a model of QDs [l4| resembling those from self- 
assembled growth 15] along crystal axis [001]. The lat- 
eral potential V(p) contains a part V s symmetric for the 



inversion p— >— p, and it may also contain an inversion- 
asymmetric part V a ^(|- We take the principal axes e x y 
of the QD to be along the crystal axes [110] and [110]. To 
construct accurate states we use a large basis set of har- 
monic oscillator wavefunctions. The lateral sizes a^ay 
are given by the curvature at the potential minimum |17| , 
which is determined entirely by the symmetric part V s . 
V a contains y nx (Vay), odd in x (y), and is parametrized 
by E x (E y ) [lj]. The effective lateral electric field in the 
ground state (ipi \— d XtV V a \<pi) oc E x , y and vanishes for 
lateral inversion symmetry. 

First we consider the two-electron wavefunctions with- 
out s-o coupling. They are obtained by diagonalizing 
hamiltonian Hq that contains the Coulomb interaction 
Uc=-^— is the dielectric constant) with a band-mixing 
correction 7 c 6(ty) 01 an d the QD potential V(r), in the 
basis of harmonic oscillator wavefunctions. These basis 
functions separate into the symmetric and antisymmetric 
sets {S 1 ,^}, {T,n^} by their permutation symmetry [l8| . 
In general, each eigenstate of Hq can be written in terms 
of functions having definite s, x, y, d symmetry, e.g.: 

T x = Tf + E X T( + E x E y Tf + E y T? , 

T 2 = T 2 + E y T 2 + E x E y T 2 + E X T 2 , (3) 

5*2 = £2 + E X S 2 + E x E y S\ + E y S 2 , 

S3 = S% + E y S^ + E x E y S§ + E X S^ , 

where T\ (T 2 ) labels the lowest triplet with a larger x (y) 
component, and S 2 (S3) labels the lowest singlet with 
a larger x (y) component. The lowest lying states are 
shown in Fig^a), where E y —Q. Higher lying states not 
shown are T3 (T4), which are the lowest d (s) symmetry 
triplets, and Si (the lowest d-symmetry singlet). The 
isotropic part of the exchange for a triplet Ti (with en- 
ergy e') and a singlet Sj (with energy ep is given by 

Jij—2(e t i —e s ^)/h 2 . Here we have chosen the energy split- 
ting between the electron ground and excited states to be 
in the range 20—45 meV; this gives an exchange splitting 
(J13 between Ti and S3) of the order 5—10 meV, in the 
range of the experiments. 

Next, the triplet-singlet mixing is generated by adding 
the s-o terms h v , h c , h B to Hq. These give a hamil- 
tonian composed of a spin-symmetric part H s that con- 
serves the total spin, and a spin-antisymmetric part H a : 

H s = H + - (hi + h 2 + ls d Pr U c xpjr)'S, (4) 

H a = - (hi -h 2 + 2j s d Pr U c x p^) ' (*i ~ « 2 ) , 

where h k =hl+h B , p,=pi-p 2 , and p c =(pi+p 2 )/2. H a 
can be written as: 

H a = X)/3« • (si - a 2 )\Ti}{Sj\ + h.c. , (5) 

fa = {Ti\hi + ls d Pr U c x p^\Sj) , 

where (3ij gives the asymmetric exchange. States of dif- 
ferent total spin \S\ are coupled via the operator s±— s 2 , 
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which is equivalent to the Dzyaloshinskii-Morya form 
jr-(siXs 2 ) [I3 ■ The asymmetric exchange can be writ- 
ten: 

/3 • (si - s 2 ) = P z (si - «|) + P 1 - ■ - si) . (6) 

The longitudinal component (3 Z conserves the total spin 
projection S z , i.e., it mixes singlets with triplets that 
have S z =0 ("longitudinal mixing"). This is equivalent to 
a precession of the total spin around axis e z (AS z =0). 
The transverse components (3 mix states with differ- 
ent total-spin projection (|AS' z |^0), equivalent with a 
total-spin precession around in-plane axis ("transverse 
mixing" ) . 
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FIG. 1: (Color online) QDs with a x —a y —a. (a) Two-electron 
energy levels in QDs with a=10 nm with one plane of symme- 
try along e x (E y =0) vs the lateral asymmetry parameter E x . 
(b) The asymmetric exchange |/xj 3 1 from Eq.® and its com- 
ponents |/ii3 V | (from s-o coupling), l/x^ ] (from spin-relative 
orbit coupling) vs the QD size a. (c) | /J.f 3 1 vs E x ^0 (E y =0) 
for several QDs. 

It is convenient to group the operators from the matrix 
element giving /Jy in Eq. (J5J into an axial vector operator 
A=Ae z and two polar vector operators P=Pe z , R=R ± : 

A = 2 7 J (d Pl V s x B P1 ) - 2 7s (d Pr U c x d Pc ) , 

P = 2 7 r (d Pl V a x d Pl ) (7) 

R = -2-{ V (e z x d p ^j + 2j B (eJ Xl - e y d v ^ . 

A and P include the vertical magnetic field from the 2D 
motion in the nearly degenerate excited states, and they 
generate (3 Z in Eq.©. R arises from the Rashba and 
Dresselhaus terms, and it generates (3- 1 . 
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TABLE I: The part of the longitudinal coupling /3y determined by A [Eq.Q]. The matrix elements Aff=(Ti"\A\S?) are 
between wavefunction components of definite symmetries [Eq.©]. 





S\ (~ s-symmctry) 


S2 (~ ^-symmetry) 


•S3 (~ y-symmctry) 


S4 (~ (/-symmetry) 


\~ s ) 


K F, \A xy 4-A yx -\- A sd -\-A ds 1 


K \ A sd -\- A yx -\-F, 2 ( A xy A- A ds W 

&y |/*42 1 ^42 tfi i ^ ■ 42 ' ^42 J J 


[^43 +-^43 1^43 1 ^43 !\ 


A| 4 + 

a x A ii T-^H 1 . 44 +- c 'x^ 1 44^ 


Ti 
(~ x) 




BxP^+A^+A^+AS] 


AH + 

ElAli+El{Aii+ElAll) 


B^A^+A^+BjCAfl+A^)] 


T 2 
(~ V) 


^[^+^+^21+^?)] 


A 22 + 


E X Ey [A^+A^ +^-23+^23] 


b«KI+a^+e2(a|i+a;|)] 


T 3 
(« d) 


A*? + 


PxIA^+A^+P^A^+Ag)] 


E y [Ail+All+El{All+A%i)} 


E*E y [AH+All+Ali+A d 3 l] 



TABLE II: The part of the longitudinal coupling % determined by P [Eq.Q]. -F > i"' 3 = <2T | P| Sj 3 ) are between wavefunction 
components of definite symmetry [Eq.J^J]. 





Si 


s 2 


S3 


Si 


Ta 


E* [PS +P^+El(P^+P^)] + 

Ey [PS+Pll +El {P^+PtD] 


PIS +El PH+El(P^+El P d J) + 
Ex E y (P 4 S 2 « +P: 2 d +P^ +P 4 t ) 


P^+P^+P^P^ + E^tH 

Ex Ey (PH+P™ +P y j 3 +Pt 3 V ) 


EX [P S i V i+PSi d +El (P» 4 S )] + 

Ey [pi: +p v i?+E 2 y (p:: +p d ? )] 


Ti 


Pii+El Pg +E 2 y (Pf? +El P^) + 
E* E y (P?f+P°? +P*° +P?° ) 


Ex [P^+PiS +E 2 y (P$ +P V 12 )] + 
Ey [P?, d +P?Z+El (P??+P?° )] 


Ex [P?i+Pt% +*%(Pg+P& )] + 
Ey [P%+P% +El (Pi 3 +P? 3 )] 


pr/+Ej p« 4 " +pj (p x t p^ )+ 
Px p H (pr 4 s +p, s r +p d i+p^ ) 


T 2 


P»l +El P 2 t + El ( P'» +El PH ) + 

^ Ey (P^ + P^+PSl +P2! ) 


Ex [P^ +Pg+B$ (PS+P 22 d )] + 
Ey [P* d +P^. +El (P%+P% )] 


e* [p^+p^+eI (p; 3 +P 23 S )]+ 

Ey [P& +P2 3 y +E 2 +P** )] 


P-!i+E 2 x P d % +E 2 y {PiS+El psd+ 

Ex E y (p« 4 s +p 2 f +p%i +p 2 f ) 


T 3 


[P 3 d i +P£ +E 2 V (P 31 d +P^)} + 

Ey [P??+P 3 1 +El(P^+Pll)] 


P d 2+E 2 x P v 32 +E 2 y (PSi+El P°») + 
Ex E v (PS+P^+P 32 +P 32 ) 


^ + PJ?+JS^ (P3" +E 2 x Pi 3 ) + 

Ex Ey (P d ?+P*° +P 33 d +P^ ) 


Px [P^+P^+Bj (P 3 X 4 3 +f 3 T )] + 

s H [p^+p 34 d +E 2 x {p^+pg)] 



Table Ogives the matrix elements between Ti (i=l, 4), 
and (j=l,4), from the spin mixing operator A in 
Eq.Q). The states are characterized by the symmetry of 
their dominant wavefunction components, e.g. S2 ~ x- 
symmetry. Table ITT1 gives corresponding results from P. 
The terms in small boxes in Table [I] are dominant and 
are independent of lateral asymmetries. All the other 
terms in Tables I and II are non-zero only for cases of 
lateral asymmetry. The central 2x2 block highlighted is 
of interest for the dynamics of X - in the "p" shell 

The matrix elements in Tables I and II can be under- 
stood by writing the operators in the basis {Tm \ Sn^}- 

A +P=H m .n ( A mn + Pmn) e z \T^) (S^ \ + h.C. The 

matrix elements Ay 3 and P^ 13 in the tables are sums of 
the matrix elements A mn and P mn with the same symme- 
try. From Eq.Q, it is seen that A mn is nonzero only for 

\Tm^){Sn \ odd both in x and in y, thus A can produce 
longitudinal mixing /3?- between two-electron eigenstates 
Ti and Sj if one of these contains a x (s)-symmetry com- 
ponent, and the other has a y (d)-symmetry part [Table 
I ■ Pmn is non-zero only for QD asymmetries (V a ^0) and 
for \Tm^)(S n ^ I odd either only in a; or only in y. Thus, P 
contributes to the longitudinal spin mixing between 
Ti and Sj if one of them has asoid component and the 
other has a x 01 y component [Table ITT] . R can be writ- 
ten as fi=E m , n #mJTr ( n 0) )(sl 0) | + h.c Results for the 
matrix elements of R mn are not given explicitly here. 
They require QD lateral asymmetry and are nonzero for 



\Tm ) (S n I odd in one of x or y. They can give transverse 
spin mixing /3,^ of states with different z spin projection. 
The degree of triplet-singlet mixing is given by the ratio 
of the asymmetric to the symmetric exchange: 

My = Tr^M, , i*tj = K^Pv/Jij ■ (8) 

We now consider QDs with different asymmetries and 
consider the longitudinal spin mixing [i z from them. This 
mixing does not have contributions from the Dresselhaus 
and Rashba couplings. 

i. QDs with lateral inversion symmetry. For them 
E x =E y ={). Examples are shown in Fig^b) and by the 
E x =0 points in FigsQIa,c) and FigEl ln sucn QDs, the 
two-electron states Ti, Sj have well-defined symmetries. 
The spin-mixing is due only to A, on the second diago- 
nal (in small boxes) in Table [I] "Pure" states of x (y)- 
symmetry such as T\ (T2) couple only to "pure" states 
of y (x)-symmetry such as S3 (S2). The first order longi- 
tudinal spin mixing of T\ (T2) is from S3 (S2), which is 
the closest in energy. T3 (the lowest d-symmetry triplet) 
couples by A to s-symmetry singlets like Si. T4 (the 
lowest s-symmetry triplet) couples by A to <i-symmetry 
singlets such as S4. 

From Fig^b) and Fig|2 (at E x =0) we can see that 
the asymmetric exchange can be a substantial fraction of 
the symmetric exchange (up to ~ 50%). Fig^b) shows 
that the asymmetric exchange is smaller for larger QDs, 
which results from larger orbits giving smaller effective 
magnetic fields in the s-o coupling. In this case /3| 2 =— f3f 3 
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because of degeneracy. The orbital momentum L z eigen- 
states -^(S 2 ±iS 3 ) are strongly coupled to -j^(Ti±iT 2 ) 
and obey AL z =0. From Fig. EJwe see that the asymmet- 
ric exchange decreases as the degeneracy of the first two 
excited states is removed by different a x and a y . In this 
case L z is not conserved. The stronger confinement along 
6y (dx>Oy) leads to J\3>J 22 , and thus to | A*x3 1 < I M22 1 - 



' 1/^13 1(10,10) 
' 1/^221(10,10) 
' 1^131(11,9) 

ltol(ii,9) 



FIG. 2: (Color online) QDs with a x ^a y . Asymmetric ex- 
change /i z in QDs with a^ll nm and a y —9 nm compared to 
the longitudinal coupling in QDs with a x =a y =10 nm. 
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FIG. 3: (Color online) Mixing of states with different spin 
projection S z in QDs with a x =a y and with a plane of sym- 
metry (E y =0, E x ^0): (a) for the mixing of T\ and 52. 
(b) fii 3 for the mixing of T\ and 5*3. 

ii. QDs with a single vertical plane of reflection. For 
this case E x ^0 and E y =0. This gives more non-zero ma- 
trix elements in TableslllandlTTl e.g. now T4 is mixed with 
53 as well as with S4. This case is illustrated in Fig2Ia,c) 
and in Fig [3 \i z for the lowest triplet is seen to decrease 
with increasing E x . For these cases, the terms propor- 
tional to E x and E\ in Tableland also the terms from 
P—2j v (d Pl VaxX-dp^mn x E x from TablelTflare nonzero, 
and they tend to cancel partially the larger terms in the 
boxes in Table|U For some triplet-singlet pairs, such as S3 
and T\ , the symmetric exchange becomes larger and thus 
their mixing decreases. Other singlet-triplet pairs can be 
degenerate, such as T 2 and S 2 in Fig^a) at E x w 1.5; 



then nonzero (3 22 leads to strong singlet-triplet mixing 
[FigEl ■ For this case, L z is not conserved. Triplets with 
(L z )^i±h can be coupled to singlets that have (L z )k^Ti. 

Hi. QDs with no vertical plane of reflection. For this 
case E X ^Q and E y ^0. Then all states in Tables [I] and 
UTI are mixed, and the degree of longitudinal spin mixing 
can be larger than in the previous cases. 

In addition to the longitudinal spin-mixing described 
above, there is also mixing that changes the spin projec- 
tion S z (transverse mixing n ± ). This arises exclusively 
from the Dresselhaus and Rashba couplings, which give 
R in Eq.JTJ). For QDs with lateral inversion symmetry, R 
mixes states which typically differ by the single-particle 
energy splitting, e.g. T\ with S\ and S4 etc.. For them 
the mixing from R is small, due to large J\\ and J14. For 
QDs with only one vertical plane of reflection, R mixes 
T\ with 5*2 or S3, which are closer in energy and there- 
fore give larger mixing. We show in FigEl this transverse 
spin- mixing for T\ and S 2 and for T\ and S3. This mixing 
occurs only for non-zero asymmetric potential (E x ^0). 
It is generally smaller than the longitudinal spin mixing 
discussed earlier, but it can become appreciable for large 
asymmetries, and it is larger in smaller QDs. 

In recent experiments on light polarization reversal in 
QDs after excited state pumping 0, , the separation 
between excited p x and p y -\ike states is small and cannot 
be resolved. One picture of this effect 0, 3 is that it 
involves mixing of electron triplets and singlets with si- 
multaneous spin flips of an electron and of the hole caused 
by axially-asymmetric e-h exchange. For this mechanism, 
however, the experimental results require large QD asym- 
metries |4j . The additional mixing of triplets and singlets 
by asymmetric e-e exchange here gives consistency with 
experiments using more realistic QD potentials. 

In addition, the present results for triplet-singlet mix- 
ing provide an alternate process for the light polariza- 
tion reversal. The angular momentum from the light can 
be given to the orbital motion of excited-state electrons. 
Then the triplet-singlet coupling given here can mix two- 
electron states that differ in their orbital angular momen- 
tum by 2H, leading to reversed light polarization. This is 
given for example in QDs with unequal lateral sizes even 
in the inversion-symmetric case (z) above. 
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B. V. Shanabrook, D. Gammon and M. E. Ware. This 
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